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Abstract 

We study three different random walk models on several two-dimensional lattices by 
Monte Carlo simulations. One is the usual nearest neighbor random walk. Another is 
the nearest neighbor random walk which is not allowed to backtrack. The final model is 
the smart kinetic walk. For all three of these models the distribution of the point where 
the walk exits a simply connected domain D in the plane converges weakly to harmonic 
measure on dD as the lattice spacing (i —)• 0. Let uj{Q,-]D) be harmonic measure for 
D, and let ujs{Q,-\D) be the discrete harmonic measure for one of the random walk 
models. Our definition of the random walk models is unusual in that we average over the 
orientation of the lattice with respect to the domain. We are interested in the limit of 
(cu5(0, •;L>) — a;(0, ■;D))/6. Our Monte Carlo simulations of the three models lead to the 
conjecture that this limit equals cm,lPd{z) times Lebesgue measure with respect to arc 
length along the boundary, where the function pd{z) depends on the domain, but not on 
the model or lattice, and the constant cm,l depends on the model and on the lattice, but 
not on the domain. So there is a form of universality for this first order correction. We 
also give an explicit formula for the conjectured density po- 
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1 Introduction 


Let D be a simply connected domain in the plane which contains the origin. We introduce a 
lattice with spacing 6 and consider a nearest neighbor random walk that starts at the origin. 
We are interested in the point where the walk exits the domain. In the limit that the lattice 
spacing 6 goes to zero, the distribution of this random point converges to harmonic measure 
which is the distribution of the point where a Brownian motion exits the domain. In this paper 
we are interested in the leading order term (in 6) in the difference between the random walk 
exit distribution and harmonic measure. Our simulations will indicate that this leading order 
term is proportional to 6. So we will refer to it as the hrst order correction. 

In addition to the nearest neighbor random walk, which we will refer to as the ordinary ran¬ 
dom walk, we study two other random walks for which the scaling limit of the exit distribution 
is also harmonic measure. One model is a nearest neighbor random walk that is not allowed to 
backtrack. In other words, at each step the walk picks with equal probability one of its nearest 
neighbors other than the site it just came from. The third model is the smart kinetic walk, also 
known as the inhnitely growing self-avoiding walk mu- At each time step the walk is only 
allowed to jump to a nearest neighbor that has not been visited before and is not a trapping 
site at that time. A site is a trapping site at time n if there is no nearest neighbor path from 
the site to the exterior of the domain through sites that are unoccupied at time n. See [S] and 
references therein for more details. 

The convergence of the exit distribution of the ordinary random walk to harmonic measure 
is a well known result. An exposition may be found in |^. A proof of the functional central 
limit theorem for the random walk without backtracking on the square lattice may be found 
in [2] . The smart kinetic walk on the hexagonal lattice is equivalent to the exploration process 
for critical percolation on the triangular lattice. In this case it has been proved that the 
exit distribution converges to harmonic measure [D m m HO]. For other lattices this is only a 
conjecture, but with good numerical support [5]. (The dehnition of the exploration process in 
the radial case, i.e., between a boundary point and an interior point, is explained in section 4.3 
of dU].) 

A priori we expect that there are two lattice effects in the hrst order correction. One is that 
the lattice introduces a length scale (the lattice spacing). The other is that the lattice breaks 
rotation invariance. Given a domain D we can introduce the lattice with any orientation with 
respect to the domain, and the scaling limit will still be harmonic measure. However, we expect 
that the hrst order correction will depend on the orientation of the lattice with respect to the 
domain. We attempt to remove this dependence by averaging over rotations of the orientation 
of the lattice with respect to the domain. Thus the random walk models that we study are in 
fact rotationally invariant, and we expect that the main ehect of the lattice is the introduction 
of a length scale. There is no doubt that this averaging over the orientation of the lattice 
changes the hrst order correction in a signihcant way, and so we are not actually studying the 
hrst order correction for the model in which the lattice orientation is kept hxed. Nonetheless, 
we see the study of this rotationally averaged model as a hrst step in the study of the hrst 
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order correction. 

Given a domain D we fix some canonical orientation of the lattice and then rotate the 
lattice about the origin by an angle a with respect to this canonical orientation. There is some 
ambiguity in how one defines the point where the random walk exits the domain. The most 
natural definition is to linearly interpolate between the steps of the walk so that it becomes a 
piece-wise linear curve in the plane, and we can then consider the first point where this curve 
intersects the boundary of the domain. We will refer to the distribution of this point as the 
discrete harmonic measure at 0 and denote it by ■; D). The symbol M indicates the 

model (ordinary random walk, random walk with no backtracking, or smart kinetic walk), and 
L indicates the lattice (square, hexagonal or triangular). This is a discrete measure on the 
finite set of sites where the bonds of the lattice with orientation a and spacing 5 intersect the 
boundary. Now we average over the orientation by defining 


a;f"(0,-;D) = 


r- 27 r 




2tt 




( 1 ) 


Because of this averaging over the orientation of the lattice, -^D) is a continuous mea¬ 

sure on dD. We let a;(0, D) be harmonic measure for D. We are interested in the difference 
•; D) — a;(0, •; D). We expect it to be of order 0(5). So we would like to compute the 
limit of [a;j^’^(0, ■; D) — D)]/5 as 5 —)■ 0. 

For each of the three models we carry out simulations on three lattices - square, hexagonal 
and triangular. Much to our surprise we find that up to an overall constant, the first order 
correction is the same for all three models on all three lattices. The constant depends on the 
model and the lattice, but not on the domain. The precise conjecture is as follows. 


Conjecture 1. For each model M and each lattice L there is a constant cm,l, and for each 
simply connected domain D with smooth boundary there is a function Pd{z) on the boundary 
such that 

hm^[a;j^’^(0,-;T>) - a;(0, •; T>)] = cm,l pniz) \dz\ (2) 

Here \dz\ is the measure on the boundary given by Lebesgue measure with respect to arc length. 
The convergence is convergence in distribution. The first order difference is a signed measure, 
so the function Pd{z) is both positive and negative. Its integral along the boundary is zero. 

The function Pd{z) is determined by the following eguation. For smooth functions g{z) on 
the boundary, 

[ PD{z)g{z)\dz\= [ ^{z)u{0,\dz\]D) (3) 

JdD JdD 

where f is the harmonic function in D with boundary values given by g, i.e., it is the solution 
of Laplace’s equation with boundary data g. The derivative is with respect to the inward unit 
normal to the domain. 
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Again, we should emphasize that the above conjecture is expected to be true only if we 
average over orientations of the lattice. If we do not average over the orientations, we expect 
that the hrst order correction depends on the orientation a. We have been deliberately vague 
about just how smooth the boundary of D needs to be in the conjecture since we do not have 
any basis for a precise characterization of the smoothness needed. 

Another possible dehnition of the random walk exit point from the domain is to take the 
hrst point of the discrete random walk outside the domain and then orthogonally project this 
point onto the boundary, i.e., take the point on the boundary that is closest to the hrst step 
the walk takes outside the domain. (We will only consider domains with smooth boundary, so 
when the lattice spacing is sufficiently small this point will be unique.) With this dehnition 
the limiting measure will still be harmonic measure. However, it is quite possible that the hrst 
order correction will be diherent. Surprisingly our simulations indicate that it is in fact the 
same. 

In the next section we give a heuristic derivation of the conjecture for the ordinary random 
walk. The following section discusses how we compute the function Pd{.z). Then we discuss 
the results of our simulations. The paper ends with some conclusions and open questions. 

2 Heuristic derivation of the conjecture 

In this section we give a non-rigorous derivation of the conjecture for the ordinary random walk. 
As we will see later, this derivation only works for the square and triangular lattices. For this 
derivation we use yet another dehnition of the random walk exit point. We run the random 
walk until it hrst hits a lattice site in D which has a nearest neighbor outside of D. Then we 
project this site to the closest site on the boundary of the domain. The resulting site on the 
boundary of the domain is the dehnition of the random walk exit point that we will use in this 
section. 

Let g he a function on the boundary of the simply connected domain D. Let / be the 
solution of the continuum Dirichlet problem 

= 0, z e D, 

fiz) = 9{z), zedD (4) 

We assume that dD and g are sufficiently smooth so that / has whatever smoothness is needed 
below. The conjecture says that 

hin^[/' g{z)us{0,\dz\]D)- [ g{z) u{0,\dz\] D)] = c [ ^{z) u{0,\dz\] D) (5) 

^^0 c> JdD JdD JdD on 

where us is the exit distribution for one of the random walk models and u is harmonic measure. 

Recall that us is the uniform average over a in [0, 27r] of us^a, where us^a uses a lattice at 
orientation a. Let Ds^a be the lattice sites in D for the lattice at orientation a, and let Bs^a be 
the sites in Ds^a with at least one nearest neighbor outside of D. For x G Bs^a dehne gs,a{x) to 
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be the value of g at the point on dD that is closest to x. (So the line from x to this point is 
perpendicular to the boundary.) Let fs^a solve 

^5,afs,a(^^') 0, Z G Ds,a 

= 95,a{^) ^ ^ Bs,a ( 6 ) 

where „ is the discrete Laplacian on the lattice with spacing S at orientation a. The discrete 
Laplacian is given by 

^ 5 ,ah{z) = ^ [h{y) - h{z)] (7) 

y:\y-z\=S 

where z is a lattice site and the sum is over the nearest neighbors y of The constant a 
depends on the lattice and is chosen so that for any sufficiently smooth function h, As^ah{z) 
converges to Ah{z) as 5 —)■ 0. (On the square lattice a = 1.) If we start a nearest neighbor 

random walk at 0 and run it until it hits Bs^a, then the average of gs,a with respect to the 

resulting measure on Bs,a is /5,a(0). Similarly, /(O) is the integral of g with respect to harmonic 
measure on the boundary. So (El) can be rewritten as 

We now dehne another solution to the discrete Laplace equation with different boundary 
data. Let gs^a he the function on Bs^a obtained by restricting / to Bs^a- Let fs,a solve 

As^afs,ci(yZ^ 0, Z G 

= 95,^ ^ (9) 

Since / is harmonic, if it is sufficiently smooth then a simple Taylor series argument shows 
that As^af is 0(6^) where p = 2 for the square lattice, p = 4 for the triangular lattice, and 
p = 1 for the hexagonal lattice. The discrete maximum principle can then be used to show 
that /(O) — /5,o(0) = 0(5^). (For example, see chapter 1 of |II].) So for the square and 
triangular lattices which have p > 2, it suffices to study the difference fs,a{.0) — /5,a(0). On 
the hexagonal lattice, the error in replacing /(O) by fs,a{0) is 0{6) which is as big as the hrst 
order correction we are trying to derive. So our derivation does not work on the hexagonal 
lattice. The difference /ya(0) —/5 ,q( 0) solves the discrete Laplace equation with boundary data 
9 s,a — 95,01- For 2; G i?ya, let Hs,a{S^,z) be the probability that the random walk (started at 0 ) 
hrst hits B^ ^ at ; 2 . Then 

Aa(0) - /ya(0) = ^ 77^,a(0, z)[g5,o,{z) - g5,a{z)] 

zeBs,a 

Recall that g 5 ,a{z) is the value of g at the point on the boundary of D which is closest to z. 
Since f = g on the boundary, we can think of g 5 ,a{z) as the value of / at this closest boundary 
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point. Recall also that gs,a{z) is the value of / at 2 ;. So a Taylor expansion of the difference of 
/ at these two points yields 

95,ai^) - g 5 ,a{.z) = -^(z) dist{z, dD) + 0{dist{z, dDf) (10) 

on 

Generically, as we move along the boundary the relation of the points 2 : G Bs^a fo the boundary 
varies. For example, on the square lattice some will have two nearest neighbors that are 
outside D, while other 2 ; will only have one. So the probability Hs,a{0, z) that the random walk 
hrsts hits Bs,a at 2 : will vary signihcantly as we move along Bs,a- 

We divide Bs^a into segments whose length is of order \/5. As the lattice spacing goes 
to zero, the number of lattice sites in each segment goes to inhnity while the length of the 
segment goes to zero. So is essentially constant on each segment. Let ^{S) denote its value 
on segment S. So we now have 

fsAO) - Aa(0) ^ ^ ffs.alo, z)dist(z, dD) 

= -V|^(S)i(S,a) Vft.40,2) (11) 


where the sum on S is over the segments and 


5(5, a) 


z)dist{z, dP) 


( 12 ) 


Since the boundary is smooth and the lengths of the segments are going to zero, the portion 
of the boundary corresponding to a segment S is becoming essentially linear. Note that the 
distance from 2 : to dD is of order 5. We conjecture that as 5 goes to zero, 5(5', a)/5 will converge 
to a limit that depends only on the angle of the tangent to the boundary with respect to the 
lattice orientation. Thus when we average over a, 5{S, a) may be replaced by some constant c 
times S. For all a, approximated by a;(0. S'; D). So we now have 


1 


'>271 


[f5,a{0) - Aa(0)] da 


-c6Y.^,(S)zz((I,S;D) 

s 


( 13 ) 


After dividing by <5 this converges to 

f (14) 

JdD dn 

So this non-rigorous argument suggests that with the dehnition of the random walk exit point 
we have used here, the constant in the conjecture will be negative. 

We end this section with a probabilistic interpretation of our conjecture for the first order 
correction. Let A be a subarc of the boundary. We consider the correction for the probability 
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that the random walk exits through A, i.e., the probability the random walk exits through 
A minus the probability a Brownian motion exits through A. Up to an overall constant, our 
conjecture is that to hrst order in 6 it is given by 


'dD 


S' 


z) a;(0, \dz\] D) 


(15) 


where / is the harmonic function with boundary data g which is 1 on A and 0 on A'^. (By 
A'^ we mean the boundary points that are not in A.) The normal derivative is approximately 
{f{z + 6n) — g{z))/6. Since the normal derivative is with respect to the inward normal, z + 6n 
is a point inside D that is at a distance S from the boundary. Let 'y^ be the curve that is a 
distance 6 inside D from the boundary. In other words, is the curve traced hj z + 6n a.s z 
runs over the boundary. Let 7 q be the part of where g = 0, and 7 ^ the part where g = 1- 
Then the above is approximately 


hi 


f{z)uj{0,\dz\-,D)+ / [f (z) \dz\-,D) 

hi 


(16) 


Note that f{z) is the probability that if we start a Brownian motion at 2 ;, then it exits 
through A. We denote this by P{z —)■ A). And 1 — f{z) is P{z —)■ A^). So the above is 


hi 


P{z ^ A)u{0,\dz\;D) - / P{z ^ A^) u{0,\dz\; D) 


'A 


(17) 


The hrst term is the probability that the Brownian motion gets within 5 of the complement of 
A but then eventually exits through A. The second term is the probability that the Brownian 
motion gets within 5 of A but then eventually exits through the complement of A. 


3 Computing the conjectured difference 

In [1] we proved that for a simply connected domain D with sufficiently smooth boundary there 
is a continuous function Pd{z) on the boundary such that for sufficiently smooth g, 

[ ^h)uj{0,\dz\]D)= [ g{z)pD{z)\dz\ (18) 

JdD dn Jgf) 

where / is the solution of Laplace’s equation with boundary data g. In this section we provide 
some details on how to explicitly compute the function Pd{z) so that we can compare it with 
the results of our simulations. 

In our simulations we always compute the cumulative distribution function (CDF) of the 
exit distribution rather than its density. So what we need to compute from our conjecture 
is the integral of Pd{z) over a subarc of the boundary. The hve domains that we use in our 
simulations have the property that a ray from the origin intersects the boundary in only one 
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point. So we can use the polar angle 0 of a point on the boundary to parameterize the boundary. 
The following derivation works for any domain with smooth boundary; one just has to take 6 
to be some parameterization of the boundary. Let ge{z) be the boundary data that is 1 on the 
arc corresponding to parameter values in [0, 9] and is 0 on the rest of the boundary. Let fe{z) 
be the solution of Laplace’s equation with this boundary data. We want to compute 

Fn{0) = j^^^{z)u{0,\dz\-D) (19) 

Let h{z) be the unique conformal map from D to the unit disc which takes the origin to 
the origin and takes the point on the boundary of D with polar angle 0 to the point 1 on the 
boundary of the unit disc. We use h to change the above integral into an integral over the 
unit circle. Letting 0 be the angle variable for the unit circle, a;(0, \dz\-,D) becomes d(f)/{27i) 
and ^{z) becomes ^-(e*'^)|h'(h“^(e*‘^))|. Here is the harmonic function in the disc with 
the following boundary data. Let z{6) be the point in dD with parameter 9, and let (j){9) be 
the polar angle of h{z{9)). The boundary data for /®(e*'^) is 1 for 0 = 0 to 0 = 0(6*), and 0 
elsewhere. We now have 

1 j'2-K 

A straightforward computation using the Poisson kernel for the disc gives 

_ sin(0(g)/2) _ 

dn^ ^ 27rsin(0/2)sin((0-0(0))/2) ^ ^ 


Thus 


r-27r 


Fd{9) = 


sin(0(0)/2) 


(27r)2 Jq sin(0/2) sin((0 - (j){9))/2) 


\h'{h-\F^))\d(l) 


( 22 ) 


We must evaluate this integral numerically. Note that it is only conditionally convergent because 
of the singularities at 0 = 0 and 0 = 0(6*). 


4 Simulations 

We have performed Monte Carlo simulations to compute the exit distribution for the three 
different walks (the ordinary random walk, the random walk with no backtracking and the 
smart kinetic walk) on three different lattices (square, hexagonal and triangular) in hve simply 
connected domains. In all cases the origin is contained in the domain. The domains are shown 
in hgure [U and dehned as follows. 


Di = {(x,?/) : ||(a;,|/) - (1,0)11 < 2}, 







D 2 = {{x,y) : -1 <y < 2 }, 

D 3 = T, 

Da = {{x,y) : \\{x,y) - {-3/4:,-l)\\ <5/2, y > -1}, 

D, = {{x,y):{x,y)-{l/2,-l/4)er} 

where T is the equilateral triangle with vertices (2,0), (—1, \/3) and (—1, —\/3), and T' is the 
equilateral triangle with vertices (3, 0), (—3/2, 3\/3/2) and (—3/2, —3\/3/2). We always start 
the walk at the origin. For all three domains the domain has been scaled so that the distance 
from the origin to the boundary is 1. 



Figure 1: The hve domains. The dot is the origin. 

We have done simulations with lattice spacings of 5 = 0.04, 0.02 and 0.01 for the three types 
of walks and the three different lattices in the hve domains. For each case we generated one 
billion samples. The average over the orientation a of the lattice is carried out as part of the 
Monte Carlo simulation. To generate each sample we hrst pick an orientation a uniformly from 
[0,27r]. Then we run the random walk model on a lattice at that orientation until it exits the 
domain. 

In our simulations we display the results using the cumulative distribution function (CDF) 
rather than the density. (Computing the density from a simulation requires taking a numerical 
derivative and so adds further uncertainty.) The differences that we plot are the difference 
of the CDF for the random walk model and the CDF for harmonic measure. We plot this 
difference as a function of the polar angle of the boundary point with respect to the origin. 
We will always refer to this function as the difference function or simply the difference. Our 
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conjecture for this function is S cm,l Fd{0) where Fr){6) is from the previous section. In the 
plots we normalize the polar angle 9 by dividing it by 27r. For all the simulations described 
here we dehne the exit point of the random walk by linearly interpolating it between lattice 
sites and looking for the hrst time this linear interpolation crosses the boundary. 

We hrst study how the overall magnitude of the difference depends on the lattice spacing 
5. Our conjecture says it should be proportional to 5. In hgure [2] we plot the norm of the 
difference vs. the lattice spacing for all three models on the square and triangular lattices for 
domain Di. The lines shown are hts using linear regression. The hgure supports our conjecture 
that the magnitude is proportional to 5. The analogous plots for the hexagonal lattice and the 
other domains are similar. In the plot RW stands for the ordinary random walk, RWNB stands 
for the random walk without backtracking, and SKW stands for the smart kinetic walk. 
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Figure 2: The L} norm of the diherence function from simulations as a function of 
the lattice spacing 5. The domain is Di. 

Next we test the part of the conjecture that asserts the constant cm,l only depends on the 
model and the lattice, not on the domain. To test this we compute the ratio of the norm 
of the simulation difference function to the norm of 5pd{z). The conjecture says that this 
ratio should be independent of the domain. Table [H shows this ratio for all three models on all 
three lattices for all hve domains. If the conjecture is true, then the values in the table should 
be the same in each row. Of course, they are not exactly the same since 5 is not quite zero 
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Di 

D 2 

D 3 


D 5 

RW hex: 

0.3996 

0.3842 

0.3833 

0.3885 

0.3641 

RW sq 

0.3611 

0.3650 

0.3704 

0.3655 

0.3429 

RW tri 

0.3668 

0.3551 

0.3626 

0.3656 

0.3539 

RWNB hex 

1.2108 

1.1954 

1.1841 

1.2051 

1.1458 

RWNB sq 

0.7625 

0.7825 

0.7670 

0.7616 

0.7376 

RWNB tri 

0.5794 

0.5745 

0.5651 

0.5699 

0.5500 

SKW hex 

1.1966 

1.1710 

1.1684 

1.1981 

1.1400 

SKW sq 

1.1632 

1.1614 

1.0951 

1.1316 

1.0803 

SKW tri 

1.0713 

1.0709 

1.0150 

1.0431 

1.0171 


Table 1: Each value in the table is the ratio of the norm of the simulation difference 
function to the norm of Sp£,{z). The row labels indicate the model and the lattice, 
and the column labels indicate the domain. 


and there are statistical errors from the simulations. The difference between the smallest and 
largest value in a given row ranges from less than 4% to a little more than 7%. 

Finally we compare the polar angle dependence of the difference functions from our simu¬ 
lations with the function predicted by our conjecture. To do this we rescale all the difference 
functions by dividing them by 6 cm,l- The conjecture says that for a given domain D, all nine 
of these rescaled difference functions should be equal to the function Fd{ 6 ). For the values 
of cm,l we use the average of the values in the table over the hve domains. We only consider 
the simulations with the smallest lattice spacing {6 = 0.01) for this comparison. Since there 
are three models being simulated on three lattices, for each domain we have nine difference 
functions in addition to the function from the conjecture. Rather than attempt to show all ten 
of these functions in a single plot, we plot the three difference functions for a single model on 
all three lattices in a single plot. For reasons of space we only show the results for the three 
domains 02 , and D^. The agreement between the simulations and the conjecture for the 
other two domains is at least as good as the agreement for these three. The resulting nine 
hgures are hguresISltofTTl Each hgure contains four curves, but it can be difficult to distinguish 
them since they are so close together. 

We have not shown any error bars in hgures [3] to [TT] to keep these hgures simple. In hgure 
[T3] we show a few error bars for the simulations of the three models on the square lattice for 
domain D 2 . The curve shown is the conjectured diherence function. The error bars for the 
other lattices and domains are similar in size. In this hgure the error bars have been rescaled by 
dividing them by 6 cm,l just as in the previous hgures. Because of this rescaling the error bars 
for the smart kinetic walk appear smaller than those of the random walk with no backtracking, 
which in turn appear smaller than those of the ordinary random walk. Before the rescaling the 
error bars are comparable in size. These error bars only represent the error from the Monte 
Carlo simulation, i.e., from the fact we can only generate a hnite number of samples. There 
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are also errors from the fact that the lattice spacing is nonzero. They are not included in the 
error bars in the hgure. 



Figure 3: Rescaled difference for the smart kinetic walk for domain D 2 . 


5 Conclusions 

We have studied numerically the exit distribution of the ordinary nearest neighbor random 
walk, the nearest neighbor random walk with no backtracking and the smart kinetic walk. We 
modify the usual dehnition of these models by also averaging over the orientation of the lattice 
with respect to the domain. For all three models the limiting distribution as the lattice spacing 
goes to zero is harmonic measure. Our simulations support the conjecture that for all three 
models on the three lattices studied, the difference between the random walk exit distribution 
and harmonic measure is, to hrst order in the lattice spacing 6, given by 6 cm,l pD{,z)\dz\ where 
the constant cm,l depends only on the lattice and the model, and the density function Pd{,z) 
depends only on the domain D. Thus there is a sort of universality for this hrst order correction. 

The most obvious open problem is to prove this conjecture; the ordinary random walk 
model is the natural model to consider hrst. In [1] the authors considered a random walk 
which takes place in the continuum rather than on a lattice. The steps of the walk are i.i.d. 
random variables which are uniformly distributed over a disc of radius 5. Note that this model 
is rotationally invariant. It was proved that the conjecture holds for this model and an explicit 
expression for the constant which is the analog of cm,l was given. This proof depended heavily 
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0.25 


conjecture 



Figure 4: Rescaled difference for the random walk with no backtracking for domain D 2 . 

on detailed estimates of the Green’s function for the model near the boundary of the domain. 
For the ordinary random walk the needed estimates on the Green’s function appear to be much 
harder. 

Our non-rigorous derivation of the conjecture for the ordinary random walk was based on 
the generator of this walk being the lattice Laplacian. The random walk with no backtracking 
is also Markovian if we enlarge the state space to include the direction the walk just came from. 
This suggests that it may be possible to give a non-rigorous derivation of the conjecture for this 
model using the generator. By contrast the smart kinetic walk is far from being Markovian, so 
a non-rigorous derivation for the conjecture for this model will require a new approach. 

In principle the non-rigorous derivation of the conjecture for the ordinary random walk gives 
a way to compute the constant, eq. (IT^ . It should be possible to compute this expression by 
Monte Garlo simulation and compare the result with the corresponding values in the table. 

Perhaps the most important question for future research is to determine the hrst order cor¬ 
rection when we do not average over the orientation of the lattice. If one hxes the orientation 
of the lattice, then the exit distribution is a discrete distribution, and this discreteness compli¬ 
cates the Monte Garlo computation of the hrst order correction in the convergence to harmonic 
measure. To avoid this one can average the orientation of the lattice only over a subinterval 
of [0,27r]. Preliminary simulations show that the difference function when we average over a 
subinterval is signihcantly different from the difference function when we average over the full 
interval. We would like to generalize our conjecture to this case of averaging over a subinterval. 
In particular, it would be very interesting to determine if the hrst order correction still has a 
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Figure 5: Rescaled difference for the ordinary random walk for domain D 2 . 

universal shape in the sense that the density function is the same for all three models on all 
three lattices except for an overall constant of proportionality. 

A Error in earlier versions 

The hrst two versions of this paper contained an error in the simulations of the ordinary 
random walk on the hexagonal lattice. The error has been corrected in this version. There are 
two changes in the results after this correction. The most signihcant change is in the values 
in the hrst row of table [H The corrected values are roughly the same size as the values in the 
second and third rows of the table for the random walk on the square and triangular lattices. 
The incorrect values for the hexagonal lattice were much larger. The other change is in the 
differencess for the ordinary random walk on the hexagonal lattice which are plotted in hgures 
EllHIIIl The overall magnitudes of these differences change signihcantly with the correction, 
but once these differences are rescaled, it is virtually impossible to distinguish the correct and 
incorrect plots. 

The error in the simulation for the ordinary random walk on the hexagonal lattice was as 
follows. At each step the walk should have three equally probable choices : backtrack to the 
site it just came from, turn left by 60 degrees or right by 60 degrees. The error was a sign error 
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Figure 6: Rescaled difference for the smart kinetic walk for domain D^. 

so that when the walk should have backtracked it instead took another step in the direction it 
had just taken. It is rather surprising that even with the error the rescaled differences exhibit 
the same universal behavior seen in all the other models. Our understanding of this is as 
follows. With the error, the simulation is actually simulating a kind of random walk on the 
triangular lattice. At each step the walk picks with equal probability one of three possible 
steps: go forward, turn left by 60 degrees, turn right by 60 degrees. In the scaling limit the 
exit distribution of this model will be harmonic measure. The simulations indicate that our 
conjecture holds for this model as well, with a different constant cm,l from the other models. 
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Figure 10: Rescaled difference for the random walk with no backtracking for domain D^. 
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Figure 12: Error bars for the three models on the square lattice. The domain is D 2 - 
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